MAX = 10000000
isp = [1 for i in xrange(MAX)]
p = []
for i in xrange(2, MAX):
    if isp[i]:
        p.append(i)
        for j in xrange(i + i, MAX, i):
            isp[j] = 0

print 'starting...'

num = 1
for i in xrange(1, 14):
    num *= i

print '13! is %d' % num

MAX2 = MAX ** 2

def phi(n):

